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Abstract 

In this paper we construct a new H(div)-conforming projection-based p-interpolation 
operator that assumes only H r (K) n H~ 1//2 (div, A)-regularity (r > 0) on the reference 
element (either triangle or square) K . We show that this operator is stable with respect 
to polynomial degrees and satisfies the commuting diagram property. We also establish an 
estimate for the interpolation error in the norm of the space H _1 ' 2 (div, K), which is closely 
related to the energy spaces for boundary integral formulations of time-harmonic problems 
of electromagnetics in three dimensions. 
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1 Introduction and main results 

This paper addresses the problem of H(div)-conforming interpolation of low-regular vector fields 
by high order polynomials. Corresponding p-interpolation operators are relevant for the analysis 
of high order boundary element approximations for time-harmonic problems of electromagnetics. 

Aiming at high-order finite element (FE) approximations of Maxwell's equations, Demkow- 
icz and Babuska [17] introduced and analysed two projection-based ^-interpolation operators 
satisfying the commuting diagram property (de Rham diagram). These are the # -conforming 
interpolation operator Ili : H 1+r (K) — > V P {K) and the H(curl)-conforming interpolation op- 
erator n™ rl : W(K) n H(curl, K) -> Vp Cd (K); here r > in both cases, K is the reference 
element (either triangle or square), V P (K) is the set of polynomials of degree < p on K, and 
"Pp^iK) is the H(curl)-conforming (first) Nedelec space of degree p (precise definitions of all 
involved Sobolev spaces and polynomial sets are given in Section f2.ll below). 

In 2D, the operators curl and div are isomorphic. The corresponding polynomial set iso- 
morphic to the Nedelec space "P v e (K) is the Raviart-Thomas (RT) space denoted by V^[K). 
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Therefore, the results of [17J related to the operator IIp Url can be used also in the H(div)- 
conforming settings (we will denote the corresponding H(div)-conforming projection-based in- 
terpolation operator by IIp lv ). In particular, given a vector field u E H r (K) n H(div, K) with 
r > 0, the interpolant u p = Hp" u G Vp^ T (K) is defined as the sum of three terms: 

a? = Ul + u p 2 + u p 3 , (l.i) 

where ui is a lowest order interpolant, u?, is the sum of edge interpolants, and U3 is the interior 
interpolant (a more detailed description of these interpolants is given in Section I2.5|) . As follows 
from [T7], the following diagram commutes: 



Hl+r(K) W(K)nil(div,K) ^ L 2 (K) 



n dh 

p 



n°_! (i.2) 



V P (K) ^ V^(K) ^ Vp~i(K), 

where EC : L 2 {K) — > V P (K) denotes the standard L 2 -projection onto the set of polynomials 
V P (K). 

The commuting diagram property, and the corresponding ^-interpolation error estimates, 
have immediate applications to the analysis of high-order FE discretisations of time-harmonic 
Maxwell's equations. In particular, these results are critical to prove the discrete compactness 
property, which in turn implies the convergence of FE approximations for Maxwell's equations, 
as well as for the error analysis (see [U El [20J, EJ [8]). We note that classical Nedelec or RT 
interpolation operators (see, e.g., [10J) are not suitable for these purposes, as they are not stable 
(with respect to the polynomial degree p) for low-regular fields and do not work equally well for 
triangular and parallelogram elements. 

When time-harmonic problems of electromagnetics are posed in infinite domains (e.g., out- 
side a scatterer), it is convenient to reformulate them as a boundary integral equation (on the 
surface of the scatterer). The energy spaces for such boundary integral equations (BIE) involve 
Sobolev spaces of negative order for both the vector field and its divergence (a typical example 
is the space H" 1 /2(div,r) in the case of a smooth (closed) surface T). Then, it is common to use 
the H(div)-conforming boundary elements (e.g., of RT type) to discretise these BIE. The fun- 
damental problem is that the underlying integral operator is not coercive, and the convergence 
analysis of the boundary element methods (BEM) requires a suitable regular decomposition of 
the energy space into the space of divergence-free vector fields and the complementary space, cf. 
|llj . In the case of Maxwell's source problem it is possible to use a decomposition, where the com- 
plementary space is regular enough even on non-smooth surfaces. Then, the H(div)-conforming 
p-interpolation operator of Demkowicz and Babuska is applicable for the convergence and error 
analysis of the p- and the /ip-BEM (see [Sill]). However, when considering the boundary integral 
formulation for the Maxwell eigenvalue problem, the orthogonal Hodge decomposition of the en- 
ergy space (see [HIE]) must be used to prove the discrete compactness property. In this case, 
the regularity issues on non-smooth surfaces affect the smoothness of the complementary space 
and prevent one from using the known H(div)-conforming interpolation operators. Hence, the 
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aim of this paper is to introduce and analyse a new H(div)-conforming p-interpolation operator, 
which is stable with respect to p and retains the commuting diagram property analogous to 
(jl.2|) . but assumes less regularity than IIp lv (namely, H r (K) n H _1 / 2 (div, X) -regularity with 

r > 0). This new interpolation operator will be denoted by IIp lv ' 5 . 

Given a vector field u £ H r (iT) n H -1 / 2 (div, K ) with r > 0, we define the interpolant 

u p = lip v '~ 5 u £ "Pp T (if) in a similar way as the interpolant u p = n£ iv u £ 'Pp T (if) (see ([TTj)): 

u p = ui + u p 2 + uf . (1.3) 

Here, ui and u?, are exactly the same as for the interpolant IT^u (see (j2Jj)J) and (|230jh 
respectively), whereas u| £ "P? 1 ' ^) is determined by solving the following system of equations: 

(div(u - ( Ul + u p 2 + uf)),divv)^_ 1/a(J0 = Vv £ V* T ' (K), (1.4) 
(u- ( Ul + u p + uf), curl $ 0> * = V<t>eP$(K), (1.5) 

where (■, ■) h~i/2(k) an d ('i')o,K denote the H -1 / 2 ^)- and the L 2 (ET)-inner products respec- 
tively. 

It is easy to see that our construction of the interpolation operator Ilp lv ' 2 is much in the 
spirit of [17] . However, we stress the use of the i7 _1 / 2 -inner product in (II. 4j) . which is more 
natural in the boundary element settings (this is in contrast to the L 2 -inner product employed 
in the definition of Hp lv , see 

flZZP ). Thus, the H~ l / 2 {K) -inner product has to be written in an 
appropriate explicit form. Of particular importance for our analysis is the following property of 
the ff -1//2 (iiT)-inner product (u, v )jj-i/2^y for a constant function v, it reduces to the L 2 {K)- 
inner product, i.e., 

(u,1)h-i/2 (k) = (u,1)o,k VueH-V 2 (K). (1.6) 

An inner product satisfying this property is presented in the Appendix (see Lemma IA.2I and 
Lemma IA.3|) . 

In the following three theorems we formulate the main results of the paper - the properties 

of the operator n p ' 2 (all proofs are given in Section [3] below) . 

The first theorem justifies the definition of the operator and states its continuity. 

Theorem 1.1 For r > the operator 

jjdiv -i . W ( K j n H- 1 /2( diV) K ) ^ L 2( K ) n H-l/2( diV; K) 

is well defined and bounded, with its operator norm being independent of p, i.e., there exists a 
constant C > independent of p (but depending on r) such that 



n div >-! 



c < C, (1.7) 
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where \\ ■ \\z is the operator norm in the space C(H r (K)Pifl 1 / 2 (div, K), L 2 (i<QnH 1 / 2 (div, K) ) . 



,div,-i , . , , „ , , . „div -i 



Moreover, the operator U p ' 2 preserves polynomial vector fields, i.e., Tl p ' 2 v p = v p for any 



The next theorem states the commuting diagram property analogous to (jl.2 
Theorem 1.2 For r > the following diagram commutes: 

H i+r( K) ^ H r (ET)nH- 1 /2( d iv,K) ^ H-V 2 {K) 



r div, 



i 



rp*"> 2 



w/iere n p : ^ l / 2 (K) — > V V (K) denotes the H l l 2 -projector. 



,div, 



The third theorem provides an error estimate for the interpolation operator U p 2 in the 
norm of the space H _1 / 2 (div, K), which is closely related to the energy space for the electric 
field integral equation (a boundary integral formulation of Maxwell's equations in 3D). 

Theorem 1.3 Ifu£ H r (div, K) with r > 0, then there exists a positive constant C independent 
of u and p such that 

|| U - rrf'^ u|| fi _ 1/2{diV!X) < c v -^i^ ||u||Hr (d i VlJ0 . (i.9) 

Remark 1.1 Using similar constructions it is possible to introduce a stable H (div) -conforming 
p -interpolation operator for even less regular vector fields u £ H r (K) n H s (div, K), r > 0, 
— 1 < s < — i. However, our proof of the commuting diagram property carries over to this 
operator if the H S (K) -inner product reduces to the L 2 -inner product for a constant function (cf. 
property (jl.6p ), which is an open problem. Although we note that such an operator would be of 
purely theoretical interest. 

Remark 1.2 Theorem 11.31 states the interpolation error estimate for sufficiently regular vector 
fields, for which one can also apply the operator IIp lv . For BIE of electromagnetics on open 
surfaces, the solution is less regular and belongs to H r (div,T) ; where r S (—^,0) and T is an 
open surface (see [151 Section 4.4] and [H Appendix A]). To obtain the error estimate for the 
corresponding BEM in this case, one can apply the (global) orthogonal projection P p with respect 
to the energy norm \\ ■ ||x- Then the estimate for ||u — P p u||x can be reduced to the estimate 
obtained in Theorem \1.3\ in the same way as in [12] or [3]. Note that using the projector P p locally 
does not guarantee the conformity of approximations (i.e., the continuity of normal components 
across interelement boundaries). Moreover, the projector P p does not satisfy the commuting 
diagram property in (jl.8p . and, thus, it is not suitable for such purposes as the convergence 
analysis and the proof of the discrete compactness. 
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The rest of the paper is organised as follows. Section [2] gives necessary preliminaries: we 
introduce the notation, recall definitions of functional spaces of scalar functions and vector 
fields, and collect auxiliary results. In particular, we give a more detailed description of the 
interpolation operators n*, IIp lv and summarise their properties (see §2.5j) . In Section [3] we prove 
the main theorems formulated above. Finally, in the Appendix we introduce some equivalent 
norms in the Sobolev spaces H r (K) and H r {K) (r = ±|), derive expressions for corresponding 
inner products, and establish the key property (jl.6p for the i? _1 / 2 -mner product. 



2 Preliminaries 

2.1 Functional spaces and polynomial sets 

In what follows, p > will always specify a polynomial degree and C denotes a generic positive 
constant which is independent of p and involved functions, unless stated otherwise. Furthermore, 
throughout the paper, K is either the equilateral reference triangle T = {x = (xi^xi); X2 > 
0, X2 < Xiv3, X2 < (1 — xi)\/3} or the reference square Q = (0, l) 2 . A generic edge of K will 
be denoted by £, and n denotes the outward normal unit vector to dK. 

We will use the standard definitions for the Sobolev spaces H r (£l) (r > 0) of scalar functions 
on $7, see, e.g., |21| (hereafter, Q is either the unit interval / = (0, 1) or the reference element K). 
The norms in these spaces are denoted by || • ||# r (n J )- For r S (0, 1) we will also need the Sobolev 
spaces H r (£l) which are defined by interpolation. We use the real K- method of interpolation 
(see [21]) to define 

H r (Q) = (L 2 (s2), flg(O)) p (1/2 < t < 1, < r < t). 

Here, Hg(fi) (0 < t < 1) is the completion of C§°(Q) in i?*(0) and we identify H^(Q) with 
F^O). Note that the Sobolev spaces H r (£l) also satisfy the interpolation property 

H r (n) = (l 2 (Q), H 1 ^)) (0 < r < 1) 

with equivalent norms. 

The L 2 -inner product and the corresponding L 2 -norm on Q are denoted by (•, -}o n and 
II • \\o,n, respectively. For r £ [—1,0) the Sobolev spaces and their norms are defined by duality 
withL 2 (f2) = H°(Q) = H°(Q) as pivot space: 

H r (Q) = (H- r {n))', H r (Q) = (H~ r {n))', 

ii I, \{u,v)o,n\ || || \(u,v)o,n\ , , 

11^11^(0)= sup -j-j — n , \m\Hr(n)= sup — . (2.1) 

0ytveH- r (n) W v \\H-r(si) o&eH-r(n) \\v\\ H -r(a) 

Note that the Sobolev spaces H r and H r on any edge I C dK are defined by using the definitions 
of the corresponding spaces on the interval I. 
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In the Appendix we consider some other expressions for norms in the Sobolev spaces H r (K) 
and H r {K) with r = ±±. We will prove their equivalence to the norms defined above, and we 
will also derive expressions for corresponding inner products. 

Throughout the paper, we use boldface symbols for vector fields. The spaces (or sets) of 
vector fields are denoted in boldface as well (e.g., H r (i£") = (H r (K)) 2 ), with their norms and 
inner products being defined component-wise. Similarly to the scalar case, the norm and inner 
product in h 2 (K) will be denoted by (•, -)o,K and || • \\o,K, respectively, which should not lead to 
any confusion. The standard notation will be used for differential operators V = (d/dxi, 8/8x2), 
div = V •, curl = Vx, and for the Laplace operator A = div V. 

Furthermore, we will use the following spaces 

H r (div,K) := {u G H r (K); divu G H r (K)}, r > 

and 

H r (div,if) :={uGH r (iO; divu G H r (K)}, rG [— 1,— |], 

These spaces are equipped with their graph norms || • ||H r (div,ii') an d || ■ llH r (div A")' respectively. 
For r = we drop the superscript in the above notation: H°(div, K) = H(div, K). 

Finally, we will need two sub-spaces incorporating homogeneous boundary conditions for the 

~ —1/2 

trace of the normal component on 8K. By Ho(div, K) (resp., H (div, K)) we denote the 
subspace of elements u G H(div, K) (resp., u G H -1 / 2 (div, K)) such that for all v G C°°(K) 
there holds 

(u, Vv) ,K + (div u, v) 0)K = 0. (2.2) 

We note that if u G H 1 ^ 2 (div, K), then identity (|2.2p holds for any v G H Z / 2 (K) by density. 
In particular, H ^(div, K) is a closed subspace of H~ 1//2 (div, K). 

Let us now introduce the polynomial sets we need. By V V {1) we denote the set of polyno- 
mials of degree <pon the interval /, and Vp(I) denotes the subset of V P (I) which consists of 
polynomials vanishing at the end points of /. In particular, these two sets will be used for any 
edge I C 8K. 

Further, Vp(T) denotes the set of polynomials on T of total degree < p, and V 2 X V ^{Q) is 
the set of polynomials on Q of degree < pi in x\ and of degree < P2 in X2- For p\ = P2 = p we 
denote V 2 (Q) = Pp )P (Q), and we will use the unified notation V P (K), which refers to V^iT) if 
K = T and to V 2 {Q) if K = Q. The corresponding set of polynomial (scalar) bubble functions 
on K is denoted by V%(K). 

Let us denote by Vf^iK) the RT -space of order p > 1 on the reference element K (see, e.g., 
m EH]), i.e., 



Vf{K) = {V p - l {K)f®^P p - l {K) = I J 



(P p 1 _ 1 (r)) 2 exP p 1 _ 1 (T) if K = T, 
V 2 P , P -i(Q)xV 2 P ~iAQ) if K = Q. 



The subset of "P? 1 \K) which consists of vector- valued polynomials with vanishing normal trace 
on the boundary 8K (vector bubble-functions) will be denoted by 7 > ?P~' (K). 
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2.2 Auxiliary lemmas 

First, let us formulate the following result, which will be used frequently in what follows. 

Lemma 2.1 The normal trace mapping u ^ u n defines a linear and continuous operator from 
H S {K) nH- 1+s (div,K) to H- l / 2+s {dK) for s G [0,±). 

Proof. Let us denote by 7t r the standard (scalar) trace operator with 7 tr : H l ~ s (K) — ► 
H l / 2 - s {8K) for s G [0, ±), and let 7t ^ 1 : H^ 2 ~ s (dK) -> H l ~ s {K) be a right inverse of 7tr . Let 
u G H s (i<Q n H~ 1+S (div, if). Taking an arbitrary w G H 1 / 2 ~ s (dK) we integrate by parts to 
obtain 



J (u-n)v da = J (div u) 7 tr 1 w dx + ^ u • V(7 tr 1 t;) <ix 



ax K K 



□ 



< ||divu||^_ 1+s(K) \\lt r l v\\ H is {K) + ||u|| h -(jc) l|V(7 ta ; 1 w)||H-(x) 

< c (|I u IIh s (x) + ll divu llH-i+ S (^)) Ibllnva-^aK)- 
Hence, u • n G H^+^dK) and we prove the continuity of the normal trace mapping: 

ii n \Id K ( u - n ) vdc7 \ 

ll u • n \\H- i / 2 + s (dK) = su p — n — n 

< C M|u[| H «(Jf) + ll divu llH-l+ s (X)) • 

We will also need the following ^-approximation result in 2D (see [21 Lemma 4.1]). 

Lemma 2.2 Let K be the reference triangle or square. Then there exists a family of operators 
{^p}; P = 1; 2, ■ ■ ■ , 7T P : H r (K) — > V p (K) such that for any f G H r (K), r > i/iere /toWs 

11/ - Vllfl*(iO < Gp-^ll/H^jo, < t < r. 
Moreover, ir p preserves polynomials of degree p, i.e., ir p f = f if f G V P (K). 

We use this result, in particular, to prove the next lemma, which provides an optimal error 
estimate for the ^-projector ILJ 1/2 : H~ l / 2 (K) -> V P {K). 

Lemma 2.3 Lei cf> G H r (K), r > —\. Then for any p > i/zere ZioWs 

ll^-D^b-i^*) < C (p + i)-^ 2 ^u\\ H r {K) . (2.3) 
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Proof. If p = then (|2,3p is trivial. Let p > 1. First, we assume that r > 0. Using the 
standard duality argument and the p- approximation result of Lemma |2, 2\ we estimate the error 
of the L 2 -projection 11° : L 2 (K) ^(iT) in the F~ 1 / 2 -norm: 

||0-n°0||^_ 1/2(x) < ||0-n>|| Oi K sup inf 11^-^11°^ 

„ II 03 — n^^llo a" 

< ||^-n^||o,K sup 

v3effV2(/£:)\{o} llVllffVa^) 

< C(p+l)- (1/2+r) |HI™- 

This estimate yields (I2.3P due to the minimization property of the if -1/,2 -projection. 

Now, let r E (— 1,0]. Assuming that 4> S H S (K) = H S (K) with some s E (0, |) and using 
the first part of the proof, we have 

110 - n-V2 0| | fi _ 1/a {K) <c{ p + i)-a/>+-) ||^[| H . (J0 . 

On the other hand, it is trivial that 



Therefore, we prove by interpolation that 

u-u-y^w^^ < c{p+ir(y^u\\ SriK) 

< C(r)(p+ir {1/2+r) H\\Hr {K ) V<PeH s (K). 
Hence, by density of regular functions in H r (K), we obtain (|2.3p . and the proof is finished. □ 
The following lemma states the inverse inequality for polynomials on the reference element K. 
Lemma 2.4 Let v p E V P (K). Then for any s, r E [—1, 1] with s < r there holds 

\\v\\ht{K) < Cp 2{r ~ s) \\v\\ H 8(K), 

where C is a positive constant independent of p. 



For r > 0, s = the proof is based on Schmidt's inequality and given in [18] for both types 
of reference elements (see Lemma 5.1 and its proof therein). By using interpolation arguments 
and induction, this result has been extended in [TO] to the full range of parameters s, r E [—1,1]. 
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2.3 The regularized Poincare integral operators 

In [TB], Costabel and Mcintosh studied a regularized version of the Poincare-type integral opera- 
tor acting on differential forms in M. n . They proved, in particular, that this operator is bounded 
on a wide range of functional spaces including the whole scale of Sobolev spaces H r (£l) (r 6 M) 
on a bounded Lipschitz domain f2 which is starlike with respect to an open ball. Moreover, the 
essential polynomial preserving property of the classical Poincare map is retained by its regular- 
ized version. Thus, the results of [16J have immediate applications to the analysis of high-order 
elements (see, e.g., [2"0| \5\ 16]). 

Let us formulate some results of [16J in two particular cases. Namely, we will define two 
Poincare-type integral operators: one operator acts on scalar functions, and the other one acts 
on divergence-free vector fields. In both cases the functions and vector fields are defined on the 
reference element K. Denoting by B an open ball in K, let us consider a smoothing function 

6eC°°(R 2 ), supp^cB, J 9(a) da = I, a =(01,02). 

B 

Then the first regularized Poincare-type integral operator R : C°°(K) — > (C°°(K)) 2 (i.e., the 
operator acting on scalar functions) is defined as Rip = , -R2) , where 

1 

R i (x) ■= J 6(a)(xi-a,i) j tip(a + t(x- a))dtda, i = l,2. 

B 

The second operator acting on vector fields is defined as follows: 
A : (C°°{K)) 2 -» C°°(K), 

1 1 
^4u(x) := J 0(a) ^(x2 — 02) J ui(a + £(x — a)) dt — (x\ — a%) J U2(a + t(x — a)) dt^j da, 
bo 

where u = (u%, U2). 

The following properties of the operators R and A are easy to check directly (see also |16} 
Proposition 4.2]): 

(Rl) R is a right inverse of the div operator, i.e., 

div{Rtp) = ip VtpeH r (K), r>0; 

(Al) if u is divergence- free, then A is a right inverse of the vector curl, i.e., 

curl(Au)=u V u £ H r (div0, K) = {u £ H r (K); div u = in K}, r > 0. 

The operators R and A satisfy the following continuity properties (see [TBI Corollary 3.4]): 
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(R2) the mapping R defines a bounded operator H r 1 (K) — ► H r (K) for any r > 0; 
(A2) the mapping A defines a bounded operator H r (K) — ► H T+1 (K) for any r > 0. 

Furthermore, the operators i? and A preserve polynomials: 
(R3) i? maps 7> p (liQ into Pj^CfiQ; 
(A3) A maps Pjp^) into V P {K). 

We will use the operators R and A to prove the following auxiliary lemma. 

Lemma 2.5 Let r > and s > r — 1. If u £ H r (K) and divu £ H S (K), then there exist a 
function ip £ H r+1 (K) and a vector field v £ H S+1 (K) suc/i £/ia£ 

u = curl-0 + v. (2.4) 

Moreover, 

l|v|| H »+i(x) < C||divu|| HS(K) and \\tp\\ H ^(K) < C \\u\\ nr{K y (2.5) 

Proof. The proof is exactly the same as for Lemma 2.3 in [5]. We use the operators R and A 
to define ip and v: 

v := #(divu) £ U S+1 (K), ^ := ^4(u - fl(divu)) £ tfminfaH-lRl^) = H r+i( K y 
Hence, due to properties (Rl) and (Al), the vector field u can be decomposed as in f)2.4|) : 

u = (u — i?(divu)) + i?(divu) = curl?/; + v. 



Inequalities (|2.5j) are then obtained by using the continuity properties of the Poincare-type 
operators and the boundedness of the divergence operator as a mapping H r (K) — > H r (K) for 
r > 0. (cf. [51 Lemma 2.3]). □ 

Remark 2.1 Note that u £ H r (K) implies that divu £ H r (K) for r > 0. T/iai is it is 
assumed in LemmafZIA that s > r — 1. 



2.4 Discrete Friedrichs inequalities 

In this subsection we improve the discrete Friedrichs inequalities of |17} Theorem 1]. This 
improvement has also become possible due to the properties of the regularized Poincare integral 
operators which were established in [16J and summarised in the previous subsection. 

Lemma 2.6 There exist positive constants C±, C% independent of p such that 
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(i) for any u G T > ^'°{K) satisfying (u, curl f)o,K = for all ip G Vp(K), there holds 

||u||o,K < Ci ||divu||^_i^; (2.6) 

(ii) /or any u G 'Pp(if) satisfying (u, curl (/7)o,if = /or a// 99 G V P (K), there holds 

||u|| ,k < Ci ||divu|| H -i(^). (2.7) 

Proof. Following the idea of [17} Theorem 1] the proof reduces to finding a continuous right 
inverse of the divergence operator within appropriate polynomial spaces. In particular, in or- 
der to prove the first statement of the lemma one needs to construct an operator T mapping 
Vp-i(K) := {ip G Vp^i(K); j R ipdx = 0} into Vp'°(K) and satisfying the following properties: 

div {Tip) = V V V G Pp-i(K), (2.8) 

||^lkjc<C||^||^_ 1(10 VVGTViW- (2-9) 
Then, given any u G Vf T '°(K) such that (u, curl <p)o t K = for all (p G Vp(K), we prove (|2.6[) : 

||u||o_ft- = min llu — curl(^||o_ft: < llu — (u — T(divu))||o_R- 

lam 

= ||T(divu)|| 0) K < C||divu||^_ 1(K) . 

Here, divu G V p -i{K) and the existence of (p G Vp{K) satisfying curl ip = u — T(divu) follows 
from two facts: 

u-T(divu) G P* T '°(iO 

and 

div(u — T(div u)) = div u — div u = 0. 

Let us construct the operator T satisfying (|2.8p . (|2.9p . Let ip G V p -i(K). Applying the regular- 
ized Poincare operator R we define v := Rip. Then v G 7 > ? r (if), due to property (R3) of this 
operator. Moreover, using property (Rl) and the fact that J K ipdx = we conclude that v • n 
has zero average along dK: 

j v-nd<7 = J divvdx = J div(Rip)dx = J ipdx = 0. 

dK K K K 

Hence, there exists a continuous piecewise polynomial <fi defined on dK such that <p\i G Vp{£) 
for any edge £ C dK and ^ = v • n on dK. Therefore, applying the polynomial extension result 
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of Babuska et al. pQ, we find a polynomial 4> € V P {K) such that 4>\dK = 4> an d there holds 



I curl ( 



\0,K 



< 



m(K) 



< C 



H 1 / 2 (dK)/I 



m 

da 



< C 

' H~ 1 / 2 (dK) 

= C || v • n|| H -i/ 2 ( 9iir ) 

< C(||v|| ,k + ||divv||^_ 1(K) ) 
|o,K + ||div(^)||^_ lw ) 



C 



(PU Lemma 2]) 
(g=v-n) 

(Lemma 12. II with s = 0) 
(v = RtP). 



Hence, using properties (Rl) and (R2) of the operator R, we obtain 

||curl^||o,tf <CWj. 1(Jf) . 



(2.10) 



Now we can define the desired operator T as Tip = Rip — curl eft. It is easy to check that 
T : V p -i{K) -> Vf T '°(K) and K2JJ holds. Making use of (j2T0]) and the continuity of the 
operator R : IJ- l {K) -> L 2 (K) (see (R2)), we also prove (l23]l . 

The proof of statement (ii) is analogous. In this case we can use the operator R: H -1 ^) — > 
~L 2 {K) for the desired continuous right inverse of div. Then R = T maps V p -i(K) into {K) 
and (|2.7p is derived similarly as above. □ 



2.5 Existing H 1 - and H(div)-conforming interpolation operators 

Let us briefly sketch the definitions and summarise the properties of the H 1 -conforming in- 
terpolation operator 11* and the H(div)-conforming interpolation operator IIp lv (see [T7J for 
details). 

Let g G H 1+r (K), r > 0. To define the interpolant H p g, one starts with the standard linear 
interpolation of g at the vertices of K: 

g\ £ Vi(K), gi = g at each vertex of K. 

Then, for each edge i C dK, we define a polynomial g^^t by using the projection 

92/ G Vp{£) : \\(g- gi)U- 92,e\\Hi/2(i) ~+ min - ( 2 - n ) 

Extending 52/ by zero onto the remaining part of dK (and keeping its notation), using some 
polynomial extension £ p from the boundary, and summing up over all edges we define 

92 ■= E G P PW' ( 2 - 12 ) 
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Finally, we define the polynomial bubble g% by projection in the ^ 1 -semi-norm 

(f 3 eP2(K): \g-(g 1 + gP + gP)\ H1{K) ^uim. (2.13) 
Then the interpolant g is defined as the sum 

Uig:=gi + fi + S&eV p (K). (2.14) 

Now we proceed to the H(div)-conforming interpolation operator. Given a vector field 

div u G -pR r 



u G W(K) n H(div,iT) with r > 0, the interpolant u p = Uf v u G "Pp T (-FQ is also defined as 



the sum of three terms: 

U P = U! + U2 + U3. 

Here, Ui is a lowest order interpolant defined as 

ui= J2 ( / u-nd<r)&, (2.15) 

where 0£ are the standard basis functions (associated with edges £) for "Pj^iT) such that 

A. - J 1 ° n 

^' n -\0 ondK\£. 

For any edge £ C <9-ftT one has 

J(u-ux)-nda = 0. (2.16) 

Hence, there exists a function defined on the boundary dK, such that 

dtp 

— — = (u — ui) • n, tp = at all vertices. (2-17) 
da 

Then, for each edge £, we define V>2 £ W by projection 

(iP\ e -4,ct>) iil/m =Q V<f>eV%(£) (2.18) 

(see Remark IA. II for the expression of (•, '}jji/2m)- Extending tp^ by zero from £ onto dK (and 
keeping its notation), we denote by ^ ^p(-^0 a polynomial extension of f/^ from (9X onto 
K, i.e., 

^eW, ^> = ^> V'lplaxv = 0. (2.19) 

Then we set 

Ur| = u^, where ^ = curl^p- (2.20) 
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The interior interpolant iig is a vector bubble function living in 7p (K) and satisfying the 
following system of equations: 

(div(u-( Ul +uf + uf)),divv) o ,^ = Vve? p RT '°(K), (2.21) 

(u - (ui + u£ + fig), curl 0)o )JC = O V0eP°(K). (2.22) 

These interpolation operators satisfy the following properties. 

Proposition 2.1 (cf. [17\ Propositions 1-3]). 

1°. For r > the operators il* : H 1+r (iT) -> fl^-fif) and H"£ iv : H r (K) n H(div,*Q -> 
H(div, K) are well defined and bounded, with corresponding operator norms independent 
of the polynomial degree p. 

2°. The operators Hp and Hp" preserve scalar polynomials in V P (K) and polynomial vector 
fields in "P^ 1 (K), respectively. 

3°. For r > 0, the diagram in \1.2\i commutes. 

The next proposition gives optimal interpolation error estimates for the operators IL^ and 
Ilp lv . These estimates are proved in [5] (see Theorems 4.1 and 4.2 therein). 

Proposition 2.2 (i) Let g € H +r (K), r > 0. Then there exists a positive constant C 
independent of p and g such that 

If - ^)>9\m{K) < c P~ r \\9\\m+r(K)- 

(ii) Let u 6 H r (div, K), r > 0. Then there exists a positive constant C independent of p and 
u such that 

|| U - Il p lv u|| H (div,i<r) < Cp~ r ||u|| H '-(div,_ft')- 

3 Proofs of theorems 

In this section we prove the main results of the paper. 

Proof of Theorem 11.11 Let u G W{K) n H -1 / 2 (div, K), r > 0. We will study each term on 
the right-hand side of (|1.3|) . Throughout the proof we denote by s a small parameter such that 
< s < min{^,r} for given r > 0. 

Step 1. Fixing an edge i C dK and using a function 
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as a test function, we integrate by parts to obtain 

J u • n da = J (u • n) cp£ da = J (div u) (j)£ dx + J u • V(j)£ dx 
l dK K K 

< ||divu||^_ 1+s ^ \\<Pt\\m- a (K) + II u IIh s (a-) I|V^||h-*(k") 

< C((f>t, s) M|u|| H r (X ) + 1 1 div u || fl"-V 2 (.K") 



Note that if divu G H~ l+S (K) then an extension to divu G H~ 1+S (K) exists but is not unique. 
However, by assumption divu G H~ l / 2 (K) C H~ 1+S {K), which is a unique extension (see [22] 
for details). Thus, ui in (|2.15|) is well defined. Moreover, since ui is a lowest order interpolant, 
we find by the equivalence of norms in finite-dimensional spaces that 

l|ui||H(div,JO - C Yl | / u - nda ^ C (|| u IIh<-(#) + ||divu||_g-_i 

Hence, due to the finite dimensionality of ui, we obtain by using Lemma |2. II 

||(u - m) • n\\ H - 1/2+sidK) < C (J|u - Ui||h«(jo + ll div ( u - u i)\\h-^+«{K) / 

< C ^||u|| H s(Ar) + [|divu[|^_i +a ^ + ||ui||n(div,iC) y 

< C (\\u\\ur{K) + ll divu ll^-i/2(^)) • (3.23) 
Step 2. From the construction of ui and from the result of Step 1 we conclude that 

(u- ui) • n € H- 1/2+s (dK), J (u- m) • nda = 0. 

dK 

Therefore, due to the isomorphism (see |17| Lemma 2]) 

^- : H x / 2 + s {dK)/W - H~ 1/2+s (dK) = {<j> G H^ 1 ^ 2+s (dK); (u, l) 0jdK = 0}, 
da 

the function V in (|2T71) is well defined, tp G H^+^dK), tp\ e G H l l 2 (t) for any edge I C dK, 
and 

X] M/llfla/a^ < C X] II^MI H V2+, W < C'||V'|| H i/2 +s(9K) < CIKu-Ui) -nll^-i/a+s^). 

(3.24) 

Hence, (I2.18P is uniquely solvable and 

Ui\\H^)<CU\t\\H^ { iy (3-25) 
15 



Furthermore, applying the polynomial extension result of Babuska et al. pQ, we find the desired 
polynomial V>2,p € V P (K) (see (|2.19p ) satisfying 

U2, P \\HHK)<C\\4hi/2 (ey (3.26) 
Thus, uf in (f2T20]) is well defined. Putting together |3T2^|) - (f3T26]) we find 

||u^||o,a- <C ^ ||curl^2, p ||o,x < C ^2 \\^2, P \\m(K) < C || (u - Ul) • n\\ H - 1/2+ s {dK y 
£CdK iCdK 

Hence, making use of (|3.23p . we obtain 

H u 2llH(div,K) = \\*%\\0,K < C (||u|| H r W + ||divil|| jj-i/a^) • (3.27) 

Step 3. The vector bubble function U3 is uniquely defined by (|1.4p - (|1.5p . To estimate the 
norms of u| and div u| we use the discrete Helmholtz decomposition 

u| = v p + curl (f) p , (3.28) 

where <j) p E V p {K) and v p G V p T '°(K) is such that (v p , curl ip) 0tK = for all (p G V°(K). 
From (jl.4p one has by using the result of Step 1 



||divu 3 5 ||£_i / 2 (X ) < C ||div(u - u 1 )||_g-_i /2( _ K -) < C y|divu||^-_i /2(K) + |divui| 

< C (||u|| H r(K-) + ||divu||^_i /2(if ^ . (3.29) 

Then, applying Lemma l2.6l fi) and recalling that divv p = divug, we find 

IIVpHci? < C ||divvp||^_ 1(if) < C ||divuf ||^-_i/ 2( ^. (3.30) 

Since (v p , curl <f) p )o,K = 0, we estimate the norm of curl (f> p by using (|1.5|) and by employing the 
results of the first two steps: 



1 1 curl p ||o, < ||u - ui - uf || ,if < C ^||u||h^(JC) + l|divu||^-_i /2(K) J . (3.31) 

Combining (|3.29p - (|3.31|) and applying the triangle inequality we obtain, by making use of de- 
composition (|3.28p . 



l u 3llo,K + ll divu 3lli?-i/2(ii-) < C f ||u|| H r(K) + ||divu|| # _ 1/2w ] . 



T div, 



The boundedness of the operator U p ' 2 (and inequality (|1.7j) ) now follows by putting together 
the results of the three individual steps and by applying the triangle inequality. 



16 



The polynomial-preserving property of the operator Hp"' 2 easily follows from its defini- 
tion. □ 

It is essential for the proof of Theorem 11.21 given below that the H~ l / 2 {K) -inner product 
satisfies (ll.6p (i.e., reduces to the L 2 (ET)-inner product for a constant function). As it follows 
from Lemma lA.31 in the Appendix, the H~ 1 / 2 (K)-mnei product given by (|A.12|) satisfies this 
property. 

Proof of Theorem 11.21 To prove the first part of the diagram, we consider u = curl g, 

g G H 1+r (K). Let us decompose H^g and Hp"' 2 u as in (|2.14p and (jl.3p . respectively. Then it 
follows from the definitions of these interpolation operators that ui = curlgi and = curlg^ 
(cf. [17J). Hence, divu = divui = divu^ = and it follows from (jl.4p that divu| = 0. 
Therefore, decomposing u| as in (I3.28|) and comparing (|1.5p with (|2.13j) . we conclude that 

u| = curl<7|. Thus, Hp"' 2 (curl g) = curl(ITpp). 

Let us prove the second part of the diagram. For any ip G V p -i(K) there exists v p G {K) 

such that divvp = (p. Therefore, decomposing IIp lv ' 2 u as in (jl.3p . we need to show that for all 
v p G Vf^iK) there holds 

div(u-Ilp lv ' ^u),divv p ^ = (div(u- (ui + uf)),divv p ) 5 _ 1/2(K) = 0. (3.32) 

Let us also decompose v p = Hp™' 2 v p G V^(K) as in f)l .3|) : 

v p = vi + vf + vf , div vi = const, div v p 2 = 0, vf G V^ % °{K). 
Then, recalling (|1.4p . applying Lemma lA. 31 and integrating by parts, we prove (|3.32p : 

divfu — nt"' 2 u),divv B ) 

V p /' P / H ~1/2 (K) 

= (div(u-ui-uf),const)^_ 1/2w + ^div(u-np V ' 2 u) , div ^fj &y2 

= (div(u — ui — Ug), const)o,_ft: = const J (u — ui — Ug) • nda = 0. 

dK 

For the last step we used the fact that u| • ti\qk = and then applied (|2.16p . □ 

For the proof of Theorem 11.31 we will need two auxiliary results regarding the new H(div)- 
conforming interpolation operator Hp"' 2 . These results are formulated in the next two lemmas: 
the first one concerns the normal trace of the interpolant Hp"' 2 u on the boundary dK, and 
the second one states some auxiliary error estimates for Hp"' 2 . 
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Lemma 3.1 Let u G W(K) n H _1 / 2 (div, if) toifft r > 0, and let uP = Up"' G ^^(K). 
JTzen /or am/ edge ^ C dK there holds 

|| (u - u?) • n||^_ lw < Cp- 1 / 2 ||( U _ U P) . n|| ff _ va(aJ0 . (3.33) 
Proof. If u G H r (if ) n H(div, if) with r > 0, then it was proved in [U Lemma 3.3] that 
|| (u - nf u) • n||^ 1(£) < Cp- 1 / 2 || (u - u) • n\\ H - 1/2(dK) . 

We note, however, that ITp lv u • n = Ilp lv ' 1 ^ 2 u • n = (ui + u|) • n on the boundary dK, and, 
as it follows from the proof of Theorem 11.11 above, ui and u| are in fact well defined for 
u G H r (if ) D H -1 ' 2 (div, K), r > 0. Therefore, the proof of Lemma 3.3 in [4] carries over to the 
case considered in this paper and inequality f)3.33j) is valid. □ 

Lemma 3.2 Let r > and s > maxj-^r — 1}. If u G W(K) and divu G H S (K), then 

||div(u- < V '^ u)|| # _ 1/2w < Cp-^ 8 ) ||divu|| H . (J0 (3.34) 
and for any e > i/iere ZioWs 

||u - Uf v '^ u\\ ,k < C (p~ r \\n\\^r [K) + e-^-^/^HdivullH.^) . (3.35) 

TTie positive constants C in (13. 34ft and (|3.35j) are independent of u and p. 

Proof. Estimate (|3.34|) is an immediate consequence of the commuting diagram property (|1.8|) 
and Lemma 12.31 

||divu- div(np v '~ 5 u)||^_i/2 W = ||divu-n~if (divu)||^_ 1/2(if) < Cp- {1/2+s) \\divu\\ HS(K) . 

Let us now prove f)3.35|) . For p = 1 this estimate follows trivially from Theorem ll.il Let p > 2. 
Using Lemma 12.51 we decompose u as follows: 

u = curl^ + v, ipeH r+1 (K), vGH s+1 (if). (3.36) 

Moreover, the norms of v and tp are bounded as in (|2.5p . Then, applying the interpolation 
operator LTp lv ' 1 and using its commutativity with IL^ (see (|1.8p ). we write 

np iv *"^ u = uf V ''^ (curl rj}) + Up"'^ v = curl(n^) + Ilf V ~^ v. (3.37) 
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Since IIp lv ' 2 is a bounded operator preserving polynomials (see Theorem II .ip . one has for any 
polynomial v p G (P P -i(K)) 2 C V™{K): 

||v-IIp lv ' 2 v|| ,x = ||v-Vp-Iip lv ' 2 (v- v p )||o,x 

- C wj nf .^2 (H v - v pIIh^) + l|div(v-v p )|| # _ 1/2w 
v P e(P P -i(K))' ! v 

< Ce" 1 inf J|v- Vp|| H1/2 + em , (3.38) 

where e G (0, ^) is fixed, e > is arbitrarily small, and for the last step we used Lemma 5 of 
|19| as well as the boundedness of the divergence operator to estimate 

||div(v-Vp)||£_i/2 W < ||div(v- v p )||^_ 1/2+e(A:) 

< Cr 1 ||div(v - v p )\\ H -i/2+e( K) < Ce~ l ||v - v p \\ n i/2+e {K y 

Applying now Lemma [2.21 componentwise and using the first inequality in (12. 5|) . we deduce from 
(13381) that 

||v - < v '"3 v|| 0) k < Ce- 1 (p - l)-( s +V2- e ) || V || HS+1 (Jf) < Ce^p-^^HdivullH.^. 

(3.39) 

On the other hand, applying Proposition I2.2f i) and the second inequality in (|2.5p we obtain 
||curl^ - ^)h,K = \i>~ ^lMm(K) < Cp~ r \M\ H i+r (K) < Cp- r \\u\\ H r {K) . (3.40) 

Combining (I3.39|) and (I3.40p we prove (13.35H by making use of decompositions (I3.36H . (I3.37H and 
the triangle inequality. □ 

Now we are in a position to prove the main interpolation error estimate. 

Proof of Theorem 11.31 For simplicity of notation we denote uP := IL P ' 2 u G V^ T (K). Let 
us consider an auxiliary problem: find uo G H(div, K) such that 

(u-u ,v)o,^ + (div(u-u ),divv)o,^ = V v G H (div, K), (3.41) 

uo • n = u p • n on dK. 

Then, using Lemma 4.8 in [12] and applying Lemmas 13.11 and 12.11 we estimate for t = —1, — I 

||u-u ||Ht+i/2(div,/o ^ c \\( u - uP ) ■ n \\m(dK) < Cp t+1/2 ||(u-u p ) • n[| H _va(sjf) 

< Cp t+1 / 2 (\\u-vP || ,K + ||div(u-u*>)|| # _ 1/2(jn ) . (3.42) 

By the triangle inequality one has 

ll U ~ U lH-l/2(div,/n - H U ~ U 0llH-V2(div,^) + H U - uP |lH-V2(div, J ft:)- ( 3 ' 43 ) 
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For the first term on the right-hand side of (|3.43p we have by using (|3.42p with t = —1: 

||u - u ||H-i/2 (diV)ir) < Cp- 1 / 2 (||u - u p ||o,^ + ||div(u - u p )||^_ 1/2w ) . (3.44) 

Now, we consider the second term on the right-hand side of (|3.43p and prove that 

||u - u1h-i/ 2((M ) < C (V 1/2 ||u - uP|| ,K + l|div(u - uP)\\ 6 . 1/2{K) ) . (3.45) 



Denote X := H 1 ^ 2 (div, K), and let X' be the dual space of X (with li 2 (K) as pivot space). 
From |13l Section 6] we know that any w G X' can be decomposed as follows: 

w = V/ + curl<7, f€HV 2 {K)/R, g G Hq(K) n H 3 ' 2 (K), 

and 



\m/ 2 (K)/R + 1 1 9 1 1 HV 2 (K) ^ C ll w llx'- (3.46) 
Hence, recalling that uq — u p G Ho(div, K) C X, we have 

||u - u p | \^-i/2 {di K) = sup -r-r = sup M— Tj 

v ; o^wgx' ll w llx' o^wgx' ll w llx' 

-(div(u - u p ),f) 0K + (u - uP,curlg) 0K 

= sup [j — [: ■ — . (3.47) 

o/wex' ll w llx' 

Let g p G V p {K). Using (jl.5p with 4> = g p and (|3.4ip with v = curlg p , we find that 

(u - u p ,curl5p)o,_ft- = (u - u p ,curl5 p )o,if - (u - u , curlg p ) 0)K = V g p G V p (K). 

Therefore, selecting g p := U p g G V p (K) and using Proposition I2.2f i). we obtain from (|3.47p 

ll u - uP |lH-i/2(div,A-) 

||div(u - uP)\\jj- 1/2(K) \\f\\ H i/2 {K)/R + ||u - u p || ,^ \g - Ulg\ H i {K) 

— SU P ii — fi 

o^wex' ll w llx' 

||div(u - U P )||H-i/ 2(g) \\f\\ H l/2(K)/R + Cp' 1 ' 2 ||U - U p || 0i ff \\g\\ H 3/2 {K) 

— SUp j"j Tj 

o/wex' ll w llx' 

< c(||div(uo-uP)|| 5 _ 1/a(JO +p- 1 / 2 ||u -u p ||o,/f). (3.48) 

Both norms on the right-hand side of (|3.48p are estimated by applying the triangle inequality 
and inequalities (j3.42[) (with t = — 1 and t = — i, respectively): 

||div(u - u p )\\fi-i/2 {K) < ||u - u ||H-i/2 (diViE -) + ||div(u - u p )||^_ 1/2w 

< C (p- 1 ' 2 ||u - vP\\ Q>K + ||div(u - h- 1/2{k) ) (3.49) 
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and 

||u -U P || 0i ^ < ||U - U P || ,A- + ||u - Uo||H(div,A-) 

< c(||u-uf|| ^ + ||div(u-uP)||^ 1/2(x) ). (3.50) 

The desired inequality in (|3.45p then follows from (|3,48p - (|3,50p . 
Now, collecting $23$ and (|3~35|) in (ET131 . we obtain 

I|U " uP |lH-l/2(div,^) ^ C (P~ 1/2 H U " uP Ho,^ + l|div(u - U P )|| fj-x/2 {K) ) ■ 

Hence, recalling that u p = Tip"' 5 u and applying Lemma 13.21 with s = r and e = |, we arrive 
at estimate (|1.9p . □ 



A Some equivalent norms and corresponding inner products in 
the Sobolev spaces H r and H r for r = ±| 

In this appendix we consider the Sobolev spaces H r and H r on the reference element K for 
r = ±tj. We will derive expressions for norms which are equivalent to those defined in Section [2TT1 
First, let us introduce some notation. 

1°. We denote by D the polyhedron (cube or triangular prism) such that D = K x (0,1). Thus 
dD = Uf =1 ti (1 = 5 if K = T and X = 6 if K = Q). Let K = T X = {(x x , x 2 , 0); (x x ,x 2 ) G 
K}, Tj = {(xi,X2,l); (xi,x 2 ) G K}, and denote K = dD\fj. Note that K is an open 
surface. We will denote by v the outward normal unit vector to dD, and we will use the 
standard notation for the gradient V and for the Laplace operator A, both acting on scalar 
functions of three variables. 

2°. Given u G H~ l l 2 (K), we denote by uk the solution of the mixed problem: find uk G 
H X (D) such that 

Au K = in D, ^0-=uonK, u K = on dD\K. 

If u G H~ l l 2 (K), then we will use the same notation as above with K replaced by K. 

3°. Given u G H 1 / 2 (dD), we denote by u its harmonic extension, i.e., the solution of the 
Dirichlet problem: find u G ^(D) such that 

Au = in D, u = u on dD. (A.l) 

4°. Given u G il x l 2 (K), we denote by u° the extension of u by zero onto dD. Thus, u° G 
H l / 2 (8D). 
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We make use of standard definitions for the norm and the semi- norm in H 1 (D): 

\\ u \\m(D) = (|MIo,d + \ u \m(D)) 1 » \ u \h\d) = l|V«||o,z>- 

Since H l l 2 (dD) is the trace space of -ff 1 (D), the norm and the semi-norm in H l l 2 (dD) can be 
equivalently written as follows 

IMUv2 (9D) - v ig {D) \\u\\ m{D) , 

U\ aD =u 

\ u \m/*(dD) - v ™l D) \ u \hhd) = ||Vfi|| ,£>. (A.2) 

U\ dD =U 

Now we can define equivalent norms in H l l 2 (K) and H 1 / 2 ^): 

Vu 5 



\U\\hI/^ K ) - \ u lifV2(SD) - 



0,D 



(A.3) 



\\u\\ H i/2 {K) ~ inf _ \\U\\hi/ H k V (A.4) 

U\ K =u 

where || • ||#i/2(#-) is defined as in (|A.3p . because K is an open surface. 

From (|X3|) one can easily derive the expression for the corresponding H l / 2 (K)-umei product. 
In fact, applying the parallelogram law twice, integrating by parts, and recalling notations 3°, 
4°, we find (see also [17] ) 

= (^.») =(«.$ I ) ^,vek^{K). (A .5) 

\ du / 0,K \ du / 0,K 

The space H- l / 2 (K) is the dual space of H l / 2 {K). We prove the following result regarding 
an equivalent norm in 

Lemma A.l For any u E H^ 1 ^ 2 ^) there holds 

[Mlij--i/2(jr) ^ [|V«a-||o,d- (A.6) 
T/ie H^ 1 / 2 -inner product corresponding to the norm on the right-hand side of (|A.6p reads as 

(u,v) H -i/ 2 ( K ) = (u,v k )o,k = (uk,v)o,k Vn, w G H~ X I 2 {K). (A. 7) 
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Proof. Using notations 2° — 4°, we integrate by parts to obtain for any u G H l l 2 (K) and any 



(\7u K ,Vv 



du 



K 



0,0 \ dv I 0,dD 
Hence, we find from (|2,ip and (|A,3P 



du 



K 



•> ■'" ) + l^ 1 -'" ) s = :"-'-:o./v- 

<9l/ /0,A" \ 6V / 0,dD\K 



u 



H-W*{K) 



sup 



0,0 



O^eff 1 ^^) lFll#i/2(/T) 



sup 



0,0 



(A. 



o,o 



Let w := uxlif. One has w G H l l 2 (K) because = on dD\K. Moreover, w° = ukIod and, 

due to the uniqueness of the solution to the Dirichlet problem (|A.lj) , we conclude that w° = uk ■ 
Therefore, 



sup 



o,o 



Vv° 



> 



0,0 



0,0 



Vw° 



Nuk\\o,d- 



(A.9) 



0,0 



On the other hand, it is easy to see that 



sup 



0,0 



Vv° 



< l|V%||o, 



o- 



(A.10) 



0,0 



Now (|A~6j) immediately follows from fO ]| - (|A.10|l . 

Using (|A.6|) together with the parallelogram law we find 



Vu K , Vuk) Vu, v G H- 1/2 (K). 

I 0,0 



(u,v) H -i /2(K) 

Hence, integrating by parts and using notation 2°, we derive (|A,7j) . 



□ 



The following lemma states an analogous result for the space H 1 < 2 (K) which is the dual 
space of H X I 2 {K). 



Lemma A. 2 For any u G H 1 ' 2 (K) there holds 



0,0 



(A.ll) 



The H~ l I 2 -inner product corresponding to the norm on the right-hand side of (|A.11|) reads as 
(u,v) s . 1/2(K) = (u, K = ((«°)*>«) n „ V ^ wG H~ l/2 (K). (A.12) 
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Proof. Let u G H- l l 2 {K). Then u° G H- l l 2 {K) C H- X I 2 {K). Using flEIJ) and COl) we have 

K«)ttf)o,id 



sup 



sup 



o^^ei/ 1 / 2 ^) ^Wm/'^R) o^weH 1 / 2 ^) ^Wm/'^R) 

\{u,V) ,k\ 



sup 



sup 



sup 



K^,w)o,x| 



o^ t ;e// 1 /2(x) vei/ 1 /2(A-) ll^ II^i/2(k) o^v&m/ 2 {K) P?f ' - WIIh 1 / 2 ^) 



V\ K =v 



\(u,v) ,k\ n I, 
SU P — = \\ u \\h-i/2(K)- 



Hence, using (|A.6|) with u replaced by u° and with K replaced by K, we prove (lA.llj) : 



H-t/ 2 {K) — W U \\H-V 2 (K) 



0,D 



Then, applying the parallelogram law, integrating by parts, and making use of notations 2°, 4°, 
we derive (fAT2|). □ 



Remark A.l The same arguments as above can be used to find equivalent norms and corre- 
sponding inner products in the Sobolev spaces on any edge I C dK . In particular, using the 
notation analogous to 3° and 4° , we have (cf. (1A.3|) . (|A.5j) ) 



Mli?l/2(£) 



Vtf 



0,K 



du c 



1/2 w ~ \~d 



n / o,e 



Vu G H y l 2 {t\ 



dv° 
dn I o 



Vu,v G H 1/2 {i). 



The next lemma states the fact that for a constant function v in (|AT2l) the H- l / 2 (K) -inner 
product reduces to the L 2 (i^)-inner product. 

Lemma A. 3 For any u G H~ l l 2 {K) there holds 

(u, l)^-i/2 {K) = (u,l)o,K- 

Proof. We have by (TATT^ 

(u, 1)h-i/ 2 (k) = ( u > <p\k)o,K, (A.13) 

where ip{x) (x = (xi,X2,xs) G D = K x (0, 1)) solves the following mixed problem (see (|A.12p 
and notations 1°, 2°, 4°): find 92 G H l (D) such that 

Aip = 0mD, f£ = 1 on Ti = K, §fc = on T; (i = 2, . . . ,1 - 1), </? = on T z . 



It is easy to see that = 1 — £3. Then (p \x = y|x 3 =o = 1 an d the assertion follows from (|A.13p . 

□ 
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